Submitted to ApJ 

Preprint typeset using IAT^X style cmulateapj v. 10/09/06 



O 

o 



MODULATIONS IN MULTI-PERIODIC BLUE VARIABLES IN THE LMC 
J. Robert Buchler 1 , Peter R. Wood 2 and R. E. Wilson 3 

Submitted to ApJ 

ABSTRACT 

As shown by Mennickent, et al. (2003), a subset of the blue variable stars in the Large Magellanic 
Cloud exhibit brightness variability of small amplitude in the period range 2.4 to 16 days as well 
as larger amplitude variability with periods of ~140 to 600 days, with a remarkably tight relation 
between the long and the short periods. Our re-examination of these objects has led to the discovery 
of additional variability. The Fourier spectra of 11 of their 30 objects have 3 or 4 peaks above the noise 
level and a linear relation of the form f a = 2(/f, — f£) among three of the frequencies. An explanation 
of this relation requires an interplay between the binary motion and that of a third object. The two 
frequency relations together with the Fourier amplitude ratios pose a challenging modeling problem. 
Subject headings: binaries: close - stars: emission-line, Be - accretion, accretion disks - galaxies: 
individual (Magellanic Clouds) 
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1. INTRODUCTION 

A search for Be stars in the OGLE and MACHO 
databases by Mennickent, et al. (2003) uncovered 30 
stellar objects that exhibit large amplitude low frequency 
light curve modulations with frequencies in the range fjj 
= 0.001 to 0.007 d" 1 (period 140 to 1000 d), in addition 
to one or two peaks in the /& = 0.05 to 0.50 d" 1 range 
(period 2 to 20 d). Furthermore they found a tight rela- 
tion between the two peaks for all the objects, with the 
form 

fb « 35/i. (1) 

Analysis led to their suggestion that these objects are bi- 
nary stars seen under high orbital inclinations, in which 
a secondary star overflows its Roche lobe onto a hotter 
primary star with the short term periodicity being di- 
rectly related to binary nature. They then conjectured 
that the low frequency is due to an elliptical disk that 
orbits the primary star. They leave open the alternative 
of precession of a tilted disk around the primary. 

Interest in these objects extends beyond their specific 
peculiarities. One can ask the general question of why 
there are very tight binaries with dimensions similar to 
those of Main Sequence stars and periods of only a few 
days, given that protostars are too large to fit within 
such systems. Clearly some kind of orbit shrinkage is 
required. If the doubly periodic binaries owe their un- 
usual behavior to hierarchically distant third stars, as 
suggested by newly detected periodicities of this paper, 
they are likely to be interesting test objects for orbital 
evolution in multiple systems. The role of third stars 
in close binary orbital evolution has been examined over 
the past several decades - recently with a rapidly de- 
veloping literature (e.g. Mazeh & Shaham, 1979, Kise- 
lcva, Eggleton, & Mikkola, 1998, Eggleton & Kiseleva- 
Eggleton, 2001, Fabrycky k, Tremaine, 2007). Simula- 
tions show that the orbit shrinkage mechanism of Kozai 
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cycles (Kozai, 1962) with tidal friction (KCTF), leads to 
production of tight binaries and may account for virtu- 
ally all ordinary binaries with periods under 10 days. A 
third star can be effective in this regard at remarkably 
large distance and low mass. Several of the cited papers 
have excellent overviews of KCTF. This general dynam- 
ical issue arises also in the problem of orbital migration 
in exoplanet systems (e.g. Mathis & Le Poncin-Lafittc, 
2009). Unfortunately, observational probes that poten- 
tially can establish relevant statistics are beset with se- 
vere difficulties, as the distant companions may be too 
dim for detection or be spatially unresolved, while reflex 
radial velocities of the inner binaries are small, with in- 
conveniently long periods. Although major progress has 
recently been made on the observational side (Tokovinin, 
et al, 2006, Pribulla & Rucinski, 2006, D'Angelo, et al., 
2006), detection of companions and measurement of their 
orbits and masses remains very difficult and degraded by 
selection effects. Therefore any new source of informa- 
tion on possible companions to very close binaries needs 
to be exploited for insights into the general problem, in- 
cluding formation, orbital evolution, and stellar evolu- 
tion. 

When examining the MACHO database for low am- 
plitude Cepheid variables (Buchler, Wood & Soszynski 
(2009)) we came across MACHO 77.7911.26, alias OGLE 
LMC SC3 274426, in which we uncovered a relationship 
among the frequencies of the highest peaks. This object 
is one of those analyzed by Mennickent, et al. (2003). 
We therefore had a closer look at all their blue variables. 
To our surprise we noticed that the frequency relation 

fa = 2(f b -f L ), (2) 

not mentioned by Mennickent, et al. (2003), clearly ex- 
ists among 11 of their 30 objects. 

Relation 2 between a pair of high frequencies and a 
much lower frequency naturally suggests two approxi- 
mately equal periodicities and associated angular mo- 
tions, leading to a model with a third star, as discussed 
further in Section 3. However the third star character- 
istics must be consistent with observed magnitudes and 
colors, so quantification of the idea is not simple and is 
still in progress. 
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TABLE 1 
The multi-periodic blue stars 



Frequencies / in d 1 , amplitudes A in mags. 



Object 


fL 


u 


h 


Mb: A l 


A a /A L 


A b /A L 


A B /A L 


Mr: A l 


A a /A L 


Ah/Ah 


Ab/Al 


77.7911.26 


0.00166 


0.12611 


0.06469 


0.0618 


0.6363 


0.2342 


0.0798 


0.0690 


0.5514 


0.1679 


0.0577 


80.6469.95 


0.00552 


0.34850 


0.17984 


0.0438 


0.2655 


0.1746 


0.2390 


0.0461 


0.2589 


0.1798 


0.1316 


1.4174.42 


0.00439 


0.25983 


0.13427 


0.1511 


0.4197 


0.0633 


0.0066 


0.1745 


0.3385 


0.0602 


0.0151 


11.9592.22 


0.00458 


0.27956 


0.14438 


0.0406 


0.2154 


0.1927 


0.1195 


0.0592 


0.1408 


0.1086 


0.1299 


80.6468.83 


0.00702 


0.41491 


0.21451 


0.0435 


0.9404 


0.1435 


0.0522 


0.0506 


0.7367 


0.1173 


0.1119 


77.8033.140 


0.00384 


0.27843 


0.14311 


0.0933 


0.2715 


0.0981 


0.1933 


0.0991 


0.2684 


0.0894 


0.1605 


212.15675.158 


0.00589 


0.38629 


0.19903 


0.1044 


0.4655 


0.1476 


0.3281 


0.1214 


0.4682 


0.1018 


0.2657 


1.3686.53 


0.00389 


0.24923 


0.12849 


0.1426 


0.2797 


0.0282 


0.1804 


0.1680 


0.2179 


0.0543 


0.1154 


11.9477.138 


0.00486 


0.30153 


0.15570 


0.0487 


0.9897 


0.1791 


0.6932 


0.0651 


0.7421 


0.0843 


0.4615 


11.9475.96 


0.00371 


0.28704 


0.14753 


0.0612 


0.5623 




0.3980 


0.0796 


0.4755 


0.0257 


0.2953 


79.5506.139 


0.00548 


0.37722 


0.19164 


0.0986 


0.2481 




0.0049 


0.1326 


0.2907 


0.0583 


0.0112 




0.03 
0.02 
0.01 
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Fig. 1. — From top down: Amplitude Fourier spectra, followed by the spectra prcwhitened successively with the 3 frequencies fi,, f a and 
/i,. Note that f a = 2(/(, — //,). Left: Object 77.7911.26 in - the largest remaining peak at the bottom is the subharmonic f a = \ j a at 
0.06305 d -1 . Middle: Object 80.6469.95 in M s - the largest remaining peak at the bottom is again the subharmonic f s at 0.17425 d . 
Right: Object 1.4174.42 in Ms; here the subharmonic is not visible above the noise level. 



2. DATA ANALYSIS 

Using the MUFRAN software (Kollath 2008) we first 
performed a Fourier Transform over the frequency range 
to 0.99 d _1 , after extirpolation 4 of the unequally spaced 
data. The light curves were then prcwhitened with the 
dominant frequency and typically 3 harmonics, after the 

4 See Press & Rybicki (1989) for an explanation of extirpolation. 



dominant frequency had been estimated with a nonlinear 
least squares fit 5 . This procedure was repeated until no 
significant frequencies remained. Some of the data were 
also independently checked with the Phase Dispersion 
Minimization method (Stellingwerf 1978) via the routine 
PDM in IRAF to confirm the detected MUFRAN frc- 

5 Prewhitening means subtracting a Fourier fit for a specified 
frequency with some number of harmonics. 
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Fig. 2. — From top down: Amplitude Fourier spectra prewhitcned successively with the 4 frequencies /r,, f a , f s 
1.3686.53 in M R . Middle: 212.15675.158 in M s . Right: 11.9475.96 in M R . 



\fa and /;,. Left: 



quencies. 

Table 1 contains our Fourier results for the MACHO 
light curve bands M B (blue) and Mr (red). Columns 2 to 
4 show the frequencies of the dominant peaks, Jl, fa and 
fb, that satisfy Eq. 2. The Mg amplitude Al of the long 
period modulation is in column 5 followed by the relative 
amplitudes A a /AL, Ab/Ai and A s /Al, where subscript 
s hereafter refers to the subharmonic frequency f s = \f a - 
The corresponding Mr amplitudes appear in columns 9 
to 12. 

For the last two objects fb stands out above the noise 
level only in Mr. The amplitude Al is always the largest 
one in both bands. For about half the objects A a is 
greater than A s , while the reverse is true for the others. 
In two objects these amplitudes are reversed between Mr 
and M B . 

Fig. 1 illustrates Fourier spectra for objects 77.7911.26 
in M R , 80.6469.95 in M B and 1.4174.42 in M B . The top 
panel shows the Fourier spectrum while the lower panels 
display the spectra after successive prewhitcnings, with 
the highest remaining peak and its three harmonics. The 
amplitudes decrease toward the noise level after the three 
prewhitenings. The largest remaining peak is at f s for 
MACHO 80.6469.95 and 77.7911.26, whereas it is hidden 



in the noise for MACHO 1.4174.42. 

When the light curves are prewhitened with the lowest 
frequency (and three harmonics), much noise can remain 
at very low frequencies for some objects, in addition to 
the ubiquitous yearly alias at w 0.0027 d _1 . This low 
frequency structure may dominate the higher frequency 
peaks, but we ignore this junk power when labeling the 
frequency peaks in terms of amplitude. For MACHO 
80.6469.95, the yearly alias remains quite prominent at 
~ 0.0027 d _1 . When the subharmonic f s peak domi- 
nates over f a wc do the prewhitcnings successively with 
peaks /i, f a , f s and finally with fb. Of course the sec- 
ond and third prewhitenings could have been done jointly 
with f s and its harmonics. Fig. 1 presents other exam- 
ples, viz. objects MACHO 1.3686.53 in M R , MACHO 
212.15675.158 in M B and MACHO 11.9475.96 in Mr. 

In MACHO 1.3686.53, f s and f b clearly stand out af- 
ter prewhitening, although they are hardly visible in the 
original Fourier spectrum. The fb peak is quite weak 
in MACHO 212.15675.158 but pops up visibly after 2 
prewhitenings. MACHO 11.9475.96 shows a weak but 
real peak at fb after 2 prewhitenings that is somewhat 
hidden by the large yearly peak and low frequency noise. 
We therefore changed the scale on the bottom 2 panels 
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These relations are displayed schematically in Fig. 3. 
We recall the observational conditions, namely (Table 1) 
that for all objects the Fourier amplitudes satisfy Al > 
A a , that both A s and Ab are smaller than A a , and that 
the ratio A s /Ab is close to unity, but can be above or 
below. 

The combination frequency /b±/l presumably is from 
radiation of the inner binary that is reprocessed ("re- 
flected") by the third object. Note that in case A the 
motion associated with frequency /r, is retrograde, while 
in case B it is prograde. 

In all of Table l's 11 objects, Jl has the highest ampli- 
tude and appears to be reasonably sinusoidal (the har- 
monics have small amplitude). This means that most of 
the light variation is generated by a third object, separate 
from, but associated with the binary pair. Mennickent, 
et al. (2003) suggest that the third object could be an 
orbiting disk or a processing inclined disk. Other possi- 
bilities involve tides in a third star on an elliptical, long 
period orbit around the inner binary. 

We prefer to separate the pure data analysis of these 
objects that has been presented in this paper from a dis- 
cussion of the physical nature and modeling that we plan 
to address in a subsequent paper. 



Fig. 3. — Schematic Fourier amplitude spectrum. Note that 
f s = s fa- Top diagram: Observed spectrum and frequency con- 
straint. Middle and lower diagrams: Two possible sets of frequency 
assignments. 

(with the yearly peak off scale). The f a peak is the 4th 
largest if one disregards the low frequency noise. 

3. DISCUSSION 

We make the reasonable assumption, as in Mennickent, 
et al. (2003), that one component is a tight binary of 
frequency fs- We also assume that frequency /l is asso- 
ciated with a third object whose nature we leave unspeci- 
fied at this stage. Note that the light curves, prewhitened 
with the frequency Jl (which has the highest amplitude) 
and folded with f a , either suggest or are compatible with 
eclipses or ellipsoidal variations. Two sets of frequency 
assignments are possible for ft,, Case A: 

fs — > /b 
fa = 2fB 

fb = f B + h (3) 

or, Case B: 

h — ► fB 

fs = fB - fL 

fa = 2(/ fl - f L ). (4) 



4. SUMMARY 

We have found a very clear relation of the form ft = 
fh + /a/2 among the dominant observed Fourier frequen- 
cies fL, fb and f a in 11 of the 30 LMC blue variables 
discussed in Mennickent, et al. (2003). This relation is 
found in both MACHO bands, except for 2 objects where 
the fb peaks seem to be hidden in the noise in Mb- Ac- 
cordingly, strong observational constraints apply to all 
11 objects, namely (a) the frequency ratio fb/ fL ~ 35 
(or fa/fh ~ 70) found by Mennickent, et al. (2003), 
(b) the linear 3-frequency relation of this paper, (c) the 
amplitude ratios of the 4 observed peaks (Table 1), and 
(d) the absence of other significant harmonics or peaks in 
the Fourier spectra. These conditions pose an interest- 
ing and serious challenge for development of a physical 
model for these 11 objects. 
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